A STUDY 1 OF THE 
CHAIN -RELAXATION 
METHOD OF 
NETWORK ANALYSIS 



BY 

ROBERT S. WILLEY 



■Library 

Navcx] 

^iiQpoiis, JVH* 9frQKiUat ^ Schc **? 



' 















A STUDY OF THE GHA2N-JIELAXATI0K 'METHOD OF lJEB/QRK AKA LYSIS 

by 

Robert 5. Wilier 



Submitted in Partial Fulfillment of the Requirements for 
the Decree of Master of Science in Electrical Engineer- 
ing at the Moore School of Electrical Engineering, Uni- 
versity of Pennsylvania 



Library 

U. S. Naval Postgraduate School 

Annapolis, Md. 



Philadelphia, Pennsylvania 



Fay, 1949 



I a 



yj- V\X<3/- vL 

W. s-7 



v • ' 

. 

f '.%• -w\ V> ■ 



VTV-idxJ 



. - doei 



■!.D' f'‘- • >h 



L •v.*/'! 



.Cl 






4*** • 






ACK;ow r LEDGiiarr 



The helpful suggestions and constructive criticism of Dr 
John G. Brainerd are sincerely appreciated. 






Contents 



Pago 

I Introduction 1 

II Presentation of the Pethod 5 

III Illustrative y 

IV Conclusions 10 

rlntes to 

Bibliography 19 



M 


















« 

rr* 

n 










A STUDY OK TIE GIAIU-JILLAXATIOI hATHOD OF KiTWOUK ANALYSIS 



I . Introduction • 

(A) In the analysis of electrical circuits of the type commonly 
encountered in radio and communication devices, a problem of major 
importance is the choice of method by which the solution of a network 
is to be obtained. There aro in general use today two principal methods 
each having certain advantages and disadvantages for the different types 
of network. A need has long existed for a single method which would be 
universally applicable, or nearly so, and which would have the added 
advantage of reducing the effort required for solution. In a recent 
article (Bibliography l), the !t Chain-Relaxation" method has been pre- 
sented as a possible fulfillment of thin need. It is the purpose of 
this paper to illustrate the application of this new method of analysis 
as a means of comparison with the more familiar procedures. 

(B) For the past half-century or more, there has existed a grow- 
ing interest in the mathematical analysis of electrical networks. As 
tho field has grown wider and more complex, and circuits have become 
more complicated, the amount of calculation and effort needed for solu- 
tion of networks has increased many times. Today, there are in exist- 
ence several electronic and electro-mechanical devices in laboratories 
and institutions which have been expressly designed for the purpose of 
analysing the more involved circuits. Such machines can do the work 
of many individual mathematicians, who could conceivably spend hundreds 

Note: Because of time lost during the period when major effort would 

normally be applied to thesis work, resulting from a serious family 
illness, this study of tho newly-introduced chain-relaxation method 
has been accepted as a master’ o thesis, in lieu of a more comprehensive 
project which had originally been chosen. 
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of man-hours on a oin^lo problem, in a natter of minutes. It is indeed 
fortunate that such advances have been made, and such analyzing mech- 
anisms will undoubtedly continue to be improved, but they are tremend- 
ously expensive items of equipment and require highly-specialized per- 
sonnel for their operation. Unless a problem is of rather great 
importance and wide application, it is not profitable to devise a 
computing mechanism for its solution and recourse must be made to clas- 
sical methods of computation. 

These classical methods of network analysis of lumped-parameter 
circuits are usually designated as the "branch currents", "mesh-" or 
"loop-currents", and "nodal" methods. They have in common the technique 
of forming sets of simultaneous linear equations on the basis of tho 
application of certain oleaantary laws of circuit behavior. Those laws 
are: Xirchhoff's first rule, which states that the sum of the currents 

approaching and leaving a node must equal zero, Kirchhoff’s second rule, 
which states that the arura of the voltage drops around a closed path 
must be equal to zero, and Ohm’s law, modified, which states that the 
current in ary branch of a network is equal to the product of the ad- 
mittance of the branch and the net voltage across the branch. As 
normally written, these laws are: 

(a) El - 0, at a node, 

(b) if : 0, for a closed path, 

(c) ZYV = 2 
or 

£Z1 = E, for a branch, 

where the definitions to be found in Guilleaiin, Bode, etc,, (see Bibli- 
ography), for loop, node, and branch are intended. In these equations I 
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is the current in a branch, V is the potential of a node at which it 
terminates, and B ia the electromotive force of any generator in the 
branch plus the electromotive forces produced in the branch by the 
action of vacuum tubes or mutual inductances. 

The primary factor in the choice of method of attack of a problem 
of other than extreme simplicity is the number of unknowns that will 
have to be determined, l or each of the three methods mentioned, cri- 
teria have been developed (Bibliography 2, 3, 4, 7) which permit de- 
termination of the number of equations which can bo expected to de- 
velop following the application of the method to a given network. If 
the number of branches is Id, and the number of nodes is n, then in the 
branch-currents method, there arc n-1 equations resulting from the ap- 
plication of KIrchhef f ' s first law to n-1 nodes plus b equations re- 
sulting from tho application of Cht: 1 s law, modified, fox' each branch. 
Therefore b + n - 1 equations result . 

For tho mesh-current s method, tho resulting number of equations is 
b - (n - l) . and for the nodal method the number of equations is (n - l). 
(A very good section on this topic is found in Bode’s (bibliography 3) 
first chapter.) 

It is evident that the branch-current method will give tho greatest 
number of equations and it is for this reason that it lias bean almost 
entirely discarded except for very simple configurations. The mesh 
and nodal methods are those in general use today, and are the ones to 
be considered hero. 

, The nesh-currente method was formulated by t axvall late in the 
last century and is perhaps more widely known than the nodal, method, 



k 



which has achieved popularity within comparatively recent years, r-ach 
has advantages and disadvantages. Ac an example of the variation which 
can occur in the number of equations. Figure (l) of Plate I shown a 
circuit willed by loop methods might require eight ©equations (assuming 
the branches shown could net be simplified) as against three by the 
nodal method. Figure (Z) of the sane plate chows a circuit requiring 
two equations by the loop method ac against seven by the nodal method. 
In general, it nay be stated that where the number of nodes is to re- 
main fixed and bronchos added or taken away, as in an experimental 
cet-up, the nodal analysis is preferable, as the number of equations 
remains unchanged, requiring only a change in rows or columns of the 
system detenidrxmt if brandies are added, or taken away. Similarly, in 
a circuit in which portions of branches are to bo varied, keeping the 
number of meshes constant, tho mesh-current method is preferable. The 
nodal method is currently finding its widest use in analyzing vacuum 
tube circuits and in saucy cases results in considerable reduction in 
the amount of labor required for solution. 

There does exist, however, one serious drawback to tho nodal 
analysis. The present procedures cannot be applied if the network 
contains an arbitrarily' located generator of negligible internal 
impedance or having regulated terminal voltage, as it becomes im- 
possible to replace tho voltage source by an equivalent current 
source as is required by tho nodal method. The manner of replace- 
ment is clearly presented in many texts (soe, for example. Bibli- 
ography 3j 7). A second disadvantage lies in the fact that present 
methods do not permit writing nodal equations directly from a circuit 
containing mutual inductance. Preliminary calculations are required. 
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A concioe presentation of th© manner of calculation is given in 
Gardner and Barnes, "Transients in Linear Systems", (Bibliography 4) . 

It is worth mentioning at this point, that this second drawback is 
possibly on the way toward removal. Heed ( Bibliography 6) has worked 
out and presented a paper in which & straightforward system of writing 
nodal equations directly froiu the network, and also of taking into 
account the presence of voltage generates of negligible internal 
impedance, is outlined* If wo assume that Reed's method is valid, 
then the mesh and nodal methods may be considered of equal merit, and 
the choice of methods is dictated by the network characteristics. As 
before, the method should be chosen which involves the lesser amount 
of calculation, or which more quickly obtains the desired result. 

From the foregoing, it can be seen that having selected one of 
the two available methods of analysis, and assuming that all possible 
use has been made of network simplification theorems prior to commen- 
cing the analysis, a certain fixed amount of computation will remain. 
If, therefore, another method of solution were developed which con- 
sistently gave fewer unknowns in the majority of problems of the type 
encountered today, its adoption should permit a great reduction in the 
amount of effort required for solution of network problems. Such a 
method has been evolved and presented by Dr. L. Tacny-Tschiassny, in 
the Journal of the Institute of Electrical Engineers (Bibliography 1). 
II. Presentation of the Lethod. 

It is not the intention to here reproduce th© article mentioned 
above, as It is published and available, but an abstract of the salient 
features will be includ©d in order that th© illustrative examples given 
to effect a partial comparison may have meaning. Let it suffice to say 
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that the author of the referenced article was inspired by Southwell 
(Bibliography 5)# who introduced the method of “residuals" as the 
basis for his relaxation methods in stress calculations. The chain- 
relaxation method is an adaptation of these methods to use in electrical 
circuit work and provides a straightforward means of solution for roost 
networks of the form cossaonly occurring in the field of radio end com- 
munication as it is today* 

The chain-relaxation method rosy be considered as a "forcing" 
nethod, in which &o few as possible arbitrary branch currents are 
assumed and from them are derived node potentials, other branch cur- 
rents, driving currents of current generators, and e.n.f . 1 3 of voltage 
generators by the use of Kirchhoff 1 s first rule and Ohm’s law, modified. 
Because of the arbitrary assumption of initial branch currents, it be- 
comes necessary at certain nodes to add "residual" currents in order to 
nake the assumed current and potential distributions physically possible. 

A very few basic rules render the actual application of the method 
simple. If the number of generators is & and the number of arbitrary 
assumptions s, them the number of residual currents (which are the un- 
knowns) will be given by? 

( d ) £ = &-£• 

A second requirement of the method is that all branches contain- 
ing voltage sources be free of impedance. This can be accomplished 
in one of two ways - by arranging for two branches in series, one 
containing the generator impedance and the other a zero impedance 
generator, or by replacing the voltage generator with its equivalent 
current generator 1 laving an impedance branch in parallel. This is 
illustrated in Figure 3, Plato I. 
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Since it would be a remote chance that the assumed, currents were 
the actual values, the procedure is repeated s - 1 tinea with different 
arbitrary values and the results tabulated. These results will also be 
incorrect but by the principle of superposition the e.m.f.'s of genera- 
tors or their driving currents nay be nade equal to their actual values 
and the residual currents nay be made equal to zero. The values of the 
£ multiplying factors will be unknown, but can bo found by solving a 
sot cf linear equations which arc fomed properly to fulfill the re- 
quirements. 

III. The procedure is best illustrated by on example which will be 
taken step by step, but, in general, one node is selected for reference 
and an arbitrary current assigned to a branch terminating thereon. Uy 
successively applying Kirchhof f 1 s first rule and Ginn' a law, modified, 
adjacent node potentials and branch currents are derived unti 1 it be- 
comes necessary to cake an additional assumption before proceeding. 
Residual currents are assigned at any node whenever the currents meet- 
ing there do not comply with Kirchoff's first rule. 

Results are most conveniently set down in tabular form. There 
are four principal points in the procedure: 

(a) Schedule of sttps. This is formed from tiie diagram of the 
circuit and indicates the direction in which the computation is to 
proceed. A symbol is used to indicate steps of special nature, such 
as an arbitrary assumption. 

(b) The operations table. The schedule of steps is 3ct down as 
Column 1 of a table. In column 2 is indicated by mathematical for- 
mulae the nature of the computation for that step, and in £ other 
columns are placed the results of the computation for the £ sets of 
arbitrary assumptions. 
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(c) The establishment of the equations. Here the results of the 
operation table are combined in such a way as to find the values of 
the s assumed currents that will render residual currents zero and 
e.m.f.'a or driving currents of generators equal to the known values. 
This step will be clearer when following an actual example, 

(d) The computation of values. The js equations are solved for the 
necessary values of the originally assumed currents and from these val- 
ues and the operation® table, the values of tho other branch currents 
and the node potentials are computed. The principle of superposition 
is employed. The expression for any quantity contained in a column 

of the operation* a table is multiplied by the correct value of the 
multiplying factor for that column and the s products ore added. 

Considerable simplification is achieved if the initially assumed 
currents are all considered unity. Tiiat is, for the first assumption, 
consider tho current equal unity in that column, and other currents 
zero, then in second colursi, second current is ono and others are 
zero, etc. This, however, is not a requisite, as any numerical value 
could be assigned, end the correct multiplying factor would be deter- 
mined when the equations were solved. 

As an illustrative example, the simple network of Plate II is 
selected and generalized impedances used in the interest of brevity. 

In practical applications, numerical values viould be inserted at any 
point at which a saving in computation could be effected. Since the 
network ha3 four independent nodes, requiring four unknowns, but may 
be solved with only throe mesh currents, the mesh analysis is chosen. 
The solution follows on succeeding plate for one of the mesh currents, 
which is in fact also a branch current. 
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On pinto (IV), the circuit is reproduced to pemit relettering to 
show the sequence of steoa by the chain-relaxation method. Table I of 
thin plate is the operations table* Arbitrary assumptions are indica- 
ted by circles around the numeral, and steps leading to the establish- 
ment of the equations by squares. The computation will be carried 
through for on© of the assumed currents which corresponds to that 
determined by neoh analysis, and it will be shown that the two are 
equal. Mode 0 is taken as reference. The solution follows the plate 
showing circuit and the two currents are shown to be equal on follow- 
ing sheets. 

Step 1 . Assign an arbitrary value of one to the current through 
element Z^. Insert 1 in column 3, end zero in column 4 of table. 

Step 2 . The potential at 2 is the voltage rice above the refer- 
ence node. Insert lx in column 3 and 0 x in column 4. 

Step 3 . No further movement can be made without a new assumption. 
Assume I^> the current to node 0 through Z^. 

Step 4 . Tho potential of node ^ is the voltage rise above refer- 
ence or Insert 0 x ^ in column 3 and 1 x in column 4. 

Tho procedure is now apparent and will not be detailed. A few 
notes are in order, however. At step 6, wo do not knew that this 
relationship is true and might be inclined to add a residual current. 
Since only one residual current will be required, £-£-£ - 1, it 
is more advisable to let these currents stand and wait to apply the 
residual at node 7, as in this way the potential of node 7 can be 
obtained from previously computed values simply. Had we added a 
residual at node 4, thin simple relationship would not hold since 
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the expression for would not contain only the assumed and derived 
values preceding but the residual as well. In step 11, we satisfy 
Kirchhoff*s first rule by the addition of the residual 1^, and in 
step 12 we set up the final condition required for a solution. The 
equations are then formed as shown on Plate (V). In order to obtain 
a solution the residual mst equal zero and V 10 nust equal known gen- 
erator voltage. The value of in column 3 is multiplied by 3^ and 
added to the value for in column 4 multiplied by 
equated to zero. Likewise the values of V 10 in column 3 and 4 are 
multiplied by and 1^ respectively and their sun equated to the 
known £, The solution for 1^ follows and is later shown to be equal 
to the corresponding current obtained from the mesh analysis. The 
solution would be completed by solving for and computing the values 
of any branch currents or node potentials from the values of columns 
3 and 4 as previously discussed by summing the respective products. 

IV. Conclusions. 

It may appear that the saving is snail in this case. Another 
example ie shown on plates CET) end and the initial equations 

arc set up. This is a generalized version of a twin-? network used 
by the author of the original article for illustrative purposes. Hero 
it is apparent that cither the mesh or nodal methods would roquire 
four unknowns, while the chain-relaxation method requires only two, 
and therefore is much to be preferred. In general, it will be found 
that the number of unknowns is at least one less than in either the 
mesh or nodal method when the circuit is of the type illustrated. 

The addition of vacuum tubes may but does not necessarily increase 
the number of unknowns in this and the nodal method, while not 
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affecting the mesh method, but this possible disadvantage is offset 
by the fact that most vacuum tubo circuits require more mesh than 
nodal equations. 

The addition of generators may work to the disadvantage of the 
chain-relaxation method, since the number of unknowns can never be 
less than the number of generators. The other methods in certain 
networks require less unknowns as current sources or generators are 
added. 

For the types of circuits cost frequently encountered, it appears 
that the chain-relaxation method can be used to advantage. The con- 
venient tabular farm lends itself readily to computation and where it 
is desired to vary elements in the network, the effects of these varia- 
tions can bo readily dot err, lined. The form of the operation table per- 
mits the computation of only such values as are needed if the complete 
solution is not required. In any event, it should be possible to deter- 
mine by inspection of a circuit diagram prior to any computation whether 
or not the chain-relaxation procedure will reduce the required number 
of unknowns below that of the better of the usual methods. 

It is this writer’s opinion that the chain-relaxation method is 
worthy of extensive trial and further investigation. As yet, no formal 
set of rules has been made available covering such variables as the 
best sequence of steps to take or the most advantageous reference node 
to select. With practice, these points should present no great diffi- 
culty, as the optimum procedure will usually be apparent, but it would 
bo of assistance to those not familiar with the method. It is to be 
hoped that further work will be forthcoming in this field which, at 
this time, holds at the very least the promise of becoming an almost 
universally satisfactory method of network analysis. 
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